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We consider damage spreading transitions in the framework of mode-coupling theory. This 
theory describes relaxation processes in glasses in the mean-field approximation which are known to 
be characterized by the presence of an exponentially large number of meta-stable states. For systems 
evolving under identical but arbitrarily correlated noises we demonstrate that there exists a critical 
temperature To which separates two different dynamical regimes depending on whether damage 
spreads or not in the asymptotic long-time limit. This transition exists for generic noise correlations 
such that the zero damage solution is stable at high-temperatures being minimal for maximal noise 
correlations. Although this dynamical transition depends on the type of noise correlations we show 
that the asymptotic damage has the good properties of an dynamical order parameter such as: 
1) Independence on the initial damage; 2) Independence on the class of initial condition and 3) 
Stability of the transition in the presence of asymmetric interactions which violate detailed balance. 
For maximally correlated noises we suggest that damage spreading occurs due to the presence of 
a divergent number of saddle points (as well as meta-stable states) in the thermodynamic limit 
consequence of the ruggedness of the free energy landscape which characterizes the glassy state. 
These results are then compared to extensive numerical simulations of a mean-field glass model 
(the Bernasconi model) with Monte Carlo heat-bath dynamics. The freedom of choosing arbitrary 
noise correlations for Langevin dynamics makes damage spreading a interesting tool to probe the 
ruggedness of the configurational landscape. 
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I. INTRODUCTION 



The theoretical understanding of the dynamical behavior of glasses is a long outstanding problem in statistical 
physics which has recently revealed new aspects on the underlying mechanism responsible of the glass transitionBcm. 
The dynamical behavior of glasses is characterized by the fast grow of the characteristic time of relaxation processes 
in the vicinity of the glass temperature T g . This increase of the relaxation time, up to fifteen orders of magnitude in 
a relatively small range of temperatures, is usually termed as the viscosity anomaly. The first consideration of this 
anomaly, the Vogel-Tamman-Fulcher law, goes back to the twenties. However, there is still no satisfying and generally 
accepted theoretical explanation for this singular behavior. 

Currently, there are two main approaches to understand the glass transition problem. One approach (the Adam- 
Gibbs-DiMarzio theoryO) focuses on thermodynamic considerations and proposes the existence of the ideal glass 
transition T g . This is a singularity where the configurational entropy of the under-cooled liquid vanishes at T g and 
a second order phase transition characterized by a finite jump in the specific heat occurs. This scenario has been 
rediscovered in the framework of mean-field spin-glasses with one step of replica symmetry breakingo. In the mean- 
field approach, T g is the temperature where configurational entropy vanishes (the so called Kauzmann temperature) 
and also replica symmetry breaks. In this paper we will denote the glass transition both by T g and T s (in this last 
case the subindex s stands for statics). The other approach relies on mode-coupling theory and describes the glass 
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transition as a strongly non-linear dynamical effect which induces long-term memory properties in the correlation 
and response functionsa. Consequence of these effects is the existence of a dynamical singularity Td where ergodicity 
breaks and correlation functions do not decay to zero. This dynamical transition T<j is a consequence of the mean-field 
character of the mode-coupling approximation. Although these two approaches are apparently different they have in 
common their mean-field character. 

One of the most distinct features of glasses is the presence of a complex free energy landscape. The viscosity 
anomaly is a signature of activated dynamics due to the existence of a rugged free energy landscape with several 
maxima and minima separated by energy barriers and saddle points which connect them. One could think that the 
existence of this type of landscape is a necessary ingredient to find the previous scenario. Well known results on 
the number of meta-stable states in spin-glasses reveal that the interesting spin-glass behavior emerges in systems 
with an exponential number of statesu. For instance, models such as the spherical p-spin interactions spin glass 
(with p > 2), the Ising p-spin interactions spin glass (with p > 3), the Edwards- Anderson model and Potts glasses 
are characterized by an exponentially large number of meta-stable states. All these models are characterized by the 
presence of quenched disorder which facilitates, analytical treatments always at the level of mean- field theory. In the 
absence of disorder, similar results are foundocj although exact calculations for the number of meta-stable states 
turn out to be more difficult. It is largely believed that an exponentially number of meta-stable states is a necessary 
condition for the existence of replica symmetry breaking. 

There have been also recent studies of exactly solvable models which, in the absence of quenched disorder, also 
exhibit glassy behaviorU'Ll (and, in particular, activated behavior of the relaxation time). These models are character- 
ized by a small number of meta-stable states. What makes them to display glassy behavior is the presence of entropy 
barriers which leads to slow dynamics even in the absence of meta-stability. Consequently, one is tempted to conclude 
that a rugged free energy landscape with a large number of meta-stable states is not essential to find glassy behavior 
but the presence of an enormous number of flat directions in phase space. I— . 

A similar conclusion was reached in the study of the Sherrington-Kirkpatrick spherical spin-glassE2l. It was found 
that an enormous number of zero modes are responsible for the slow dynamics found in this model. Although that 
model does not have an exponential large number of meta-stable states and does not show activated dynamiGS-.it 
displays glassy behavior at low temperatures due to the existence of flat directions around the meta-stable stateald. 

A possible way to investigate the existence of flat directions in a rugged free energy landscape is the study of the 
damage spreading. Damage spreading consists in the study of the dynamical evolution of the distance D{t) (to be 
defined later) between two systems configurations evolving under the same dynamical rules and differing only in their 
initial condition (the so called initial damage). 

The study of the damage spreading problem (hereafter referred as DS) was proposed by Kaufmann-in the 60's 
for the study of the propagation of mutations in the genotype in the biological growth of individualsli3. That is, 
how a small perturbation in the genotype (microscopic level) manifests itself in the long-time term in the fenotype 
(macroscopic level). Since then, such problem has received considerable attention in the framework of statistical 
physics, particularly in the middle 80'sE2l. Almost fifteen years ago it was realized DS could be a powerful tool to 
distinguish different dynamical regimes in disordered systems, such as spin glasses£j'E£l. Variants of damage spreading 
phenomena have also been proposed-to-numerically investigate equilibrium correlation functions in generic statistical 
systemsE^I and lattice gauge theoriesED'ta. However, the initial enthusiasm and exciting perspectives in the research of 
this problem decreased in subsequent years after realizing that this transition was dependent on the type of dynamics 
used. So the existence of the DS transition could have nothing to do with the presence of a thermodynamic phase 
transition. Physicists then started to systematically investigate the DS in well known ordered systems such as the Ising 
model. In particular, much work haSrbeea devoted to the study of the one-dimensional Domany-Kinzel automatonlia 
and the one-dimensional Ising modelEj'EU. The question of the non-universality of damage spreading has been also 
emphasized in the context of non-equilibrium phenomena such as domain growth by Graham, Hernandez-Garcia and 
GrantS. 

So the question remains whether this transition has a truly physical meaning or not. In this direction, Hinrichsen 
and Domany tried to give a precise dynamic- independent definition for DS. To define a damage spreading phase one 
must consider all possible dynamical procedures which lead the system to thermal equilibrium. For discrete systems 
with small number of nearest neighbors this definition can be implemented but not in the general case (for instance, 
continuous systems) where an infinity of dynamical rules can be always implemented. 

The purpose of the present work is to present a detailed study of the DS in a model with a rugged free energy 
landscape with an exponentially large number of meta-stable states. In particular, we will study the DS in the p-spin 
spherical spin glass, an exactly solvable model for the glass transition which is described by the Adam-Gibbs-DiMarzio 
scenario and whose dynamics is described by the mode-coupling equations. To be more specific, we will study damage 
spreading for Langevin dynamics, the simplest dynamics which is continuous in time and satisfies ergodicity and 
detailed balance. It must be stressed that although there are very few works on the DS problem using this dynamics 
(StarioloES and Graham et alS), the majority of theoretical works in DS have considered discrete dynamics (in 
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discrete systems). 

We suggest that damage spreading can be used as a dynamical method to evidenciate the existence of a large number 
of flat directions or saddle points in phase space. Also the existence of a dynamical transition will be shown. Although 
we will check that damage spreading transitions are strictly non-universal still it is possible to use the asymptotic 
distance to define an order parameter for this dynamical transition. We also anticipate that by considering correlations 
between the noises of the two evolving systems an infinity of dynamical transitions can be obtained. For Langevin 
dynamics, the case of maximally correlated noises has a particularly interesting physical meaning. 

Because our study considers the DS in the mode-coupling theory for glasses it is expected to be generally valid for 
Langevin dynamics in systems with a rugged free energy landscape such as realistic glasses. Although the DS transition 
is non-universal and depends on different dynamical rules (or cross-correlations between the stochastic noises), we 
believe that this transition gives an interesting information on the free energy landscape and could be investigated 
in structural glasses. Being a signature of the existence of saddle-points in phase space (i.e. points which separate 
stable and unstable phase space directions) it is natural to expect that real glasses are good systems to manifest these 
effects—This consideration makes our results more attractive from the viewpoint of numerical simulations of realistic 
glassescj. Nevertheless we call the attention to the reader that some of our claims in this paper are not generally 
proven (such as the connections between damage spreading transitions and saddle points in phase space) and the 
present reasarch should be seen as a first step for a better understanding of some of these questions. 

The contents of the paper are divided as follows. The second section is devoted to general considerations and 
definitions about the DS problem. Section III demonstrates the existence of a dynamical transition T for DS. Section 
III is divided in three subsections. Subsection III. A describes the mode-coupling equations for the p-spin spherical spin 
glass, starting from a random initial configuration for different correlations between the noises. This subsection also 
describes the different numerical methods we have used to analyze the mode-coupling equations for different cases. 
Subsection III.B analyzes the DS problem starting from an equilibrium initial condition. Subsection III.C analyzes 
the DS in the presence of asymmetric interactions. Section IV presents an analysis of the damage spreading for glassy 
models with discrete dynamics. Concretely, we study the DS transition in the Bernasconi model with heat-bath 
dynamics. Finally section V presents the conclusions. Two appendices are devoted to some technical issues. 



II. GENERAL CONSIDERATIONS ABOUT DS 



In the most general framework the DS problem can be stated as follows. Consider a dynamical system described 
by a generic variable C which denotes a given configuration. Suppose that the system evolves under a deterministic 
dynamical rule F. For sake of simplicity we take a continuous time dynamics. The equation of motion reads, 

C t = F(C t ) . (1) 

On top of the configuration C and the dynamical rule F we also need to define a distance in the phase space of 
configurations (for instance, a Hamming distance for spin systems). This distance D needs to satisfy the usual good 
properties, in particular D(Ct,Ct) = at all times. Suppose that we take two initial configurations Co, Co with initial 
distance -Do = D(Cq,C'o) and consider the generic equal times distance, 

D(t) = D(C t) C' t ) (2) 

where C(t) and C'(t) start from configurations Co and C'o at time and evolve under the same dynamical rule F 
eq.(Q). Our main interest is to investigate the value of the asymptotic long-time distance -Doc , 

Deo = lim D{t) (3) 

t — >oo 

Note that D(t) = if Do = 0. Quite generally the asymptotic distance Doc will be a function of the type of initial 
configurations Co, C'o as well as their initial distance -Do- The dependence of Doc on those parameters is governed by 
the dynamical properties of the deterministic rule, such as chaotic properties and Lyapunov exponents. 

One could extend this general problem to stochastic systems, i.e. dynamical systems which evolve in the presence 
of a stochastic noise. Let us consider two systems described by the configuration variables C t ,C' t at time t which 
evolve following a Langevin dynamics, 

Ct = F{C t )+ Vt ; C> t = F(C' t )+r)' t (4) 

where F(C) is a external force (which can eventually derive from a potential although this is not a necessary 
condition) and rjt, rf t are external white noises uncorrelated in time with variance 2T. Suppose now we make C, C evolve 
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following eq.([|) with the same realization of the stochastic noise and starting from two different initial conditions. 
We are thinking of noises statistically identical which coincide at equal times (i.e r/ t — rf t ) when C f = C't. So we 
choose (r]i(t)rij(s)) — = 2T6(t— s)Sij and cross-correlations {r)i{t)r}' (s)) = 2TfC(Ct,C't)S(t — s)S{j where 

K is a generic function which satisfies the properties K.(C,C) = 1 and —1 < K.(C,C) < 1, VC,C. In the presence of 
stochastic noise, again D(t) = if D(0) — 0. Note that the role of the correlations is irrelevant for the evolution of 
the independent systems C,C but crucial for their correlations and the equal times distance (||). Different choices 
of the cross-correlation IC for Langevin dynamics are the analogous of different dynamical rules in discrete dynamics 
such as Monte Carlo (these rules could be Glauber, Metropolis or heat-bath among others). 

Now we are interested in the asymptotic long-time distance eq. (||) . Quite generally, will be a function of the 
type of initial condition (for instance, random or stationary), the initial distance Dq, the intensity of the noise T, and 
the cross-correlator IC. For the case in which r\ t = v( t (K = 1), we will find that there is a dynamical phase transition 
at a finite temperature To below which the asymptotic distance is different from zero. The origin of this dynamical 
transition can be explained in plain physical words. In eq.(Q) there is competition between two different terms. On 
the one hand, the force term F(Ct) propagates the error (or damage) in the initial configuration. Instead, the noise 
•qt acts in the same way in both systems smearing out possible differences in the initial condition. In other words, the 
stochastic noise is the synchronizing force which tries to make both evolving configurations to merge in time while 
the force term amplifies the initial damage playing the role of a noise. This argument only applies if fC = 1. In the 
general case — 1 < K. < 1 the noise does not necessarily synchronize both systems and its effect is similar to that of 
the force. In this case, the asymptotic distance will be also a function of the cross-correlation IC (i.e. the dynamical 
rule). 

To understand better the role of the cross correlation IC let us consider as a starting point the simple problem of a 
particle which moves inside an harmonic potential V(x) = \x 2 following a Langevin dynamics, 
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dx{t) _ dV{x) 
dt dx 



■77(f) (5) 



where rj{t) is a stochastic white noise of variance 2T. The configuration C corresponds to the position x of the 
particle and we define a distance between two configurations x, y as D(x,y) = (x — y) 2 . Take now two identical 
particles x, y and make them follow eq.(j^) both with the same stochastic noises ?y, rf and cross-correlation IC(x,y). 
For simplicity we will take a x, y symmetric cross-correlation IC — JC(D). If < .. > stands for average over dynamical 
histories then the distance D(t) =< (x(t) — y(t)) 2 > satisfies the following equation, 

¥?j& = -2D(t)+4T(l-lC(D)) . (6) 

This equation has several stationary solutions depending on IC. Obviously D = is a stationary solution (remember 
JC(0) = 1) which implies /C'(0) < 0. Only if -1/2T < fc'(0) < the solution D = is stable. On the other hand, 
for /C'(0) < — I/2T the solution D = is unstable and D(t) converges to another stationary solution (which can be 
shown it is always stable). So, there is a dynamical transition at T c = 1/(2|/C'(0)|) where the asymptotic distance 
Doo changes from zero (T < T c ) to = D* where D* satisfies the identity D* — 2T(1 — JC(D*)). The asymptotic 
distance is then given by, — limt_,oQ D(t) — D*Q(T* — T) independent of the value of the initial distance Do 
between the two particles. So when /C'(0) does not vanish, already for the simple harmonic oscillator there is more 
than one stationary solution. In other words, the effect of the cross-correlation term JC(D) manifests through the 
appearance of more than one stationary solution. As the reader can imagine the discussion turns out to be more 
difficult for other more complex potentials. 

Although this dependence of the asymptotic distance on the cross-correlation K, is an intrinsic property of 
damage spreading it does not necessarily imply that this kind of dynamical phase transitions do not give any relevant 
information on the physical properties of the system. What this really means is that the results concerning DS 
may depend on particular forms of cross-correlation function between the noises (similarly as happens for discrete 
dynamics where different dynamical rules yield different results). Nevertheless, generic results for DS may be obtained 
for correlators which satisfy quite general conditions (as we will see later). 

The major part of work in this paper will be devoted to the study of maximal cross-correlations, i.e IC = 1 although 
the results are extrapolable to more generic cross-correlators. For Langevin dynamics, the case IC = 1 is particularly 
appealing for two reasons. 

On the one hand, it followed from the simple example of the harmonic oscillator that JC = 1 implies KJ = 
everywhere. Below we will argue that this observation holds also for more complex potentials. More specific, we 
argue that in the case K, = 1 there is a single stationary solution = for any finite system and for any confining 
potential which diverges in the boundaries (i.e V(x) — > oo when x — > ±oo). Clearly, the harmonic oscillator is a 
trivial case where the asymptotic distance always goes to zero independently of the temperature T and how far is x 
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from y at t = 0. In order to justify our assertion let us take a more complex potential of two wells separated by a 
finite barrier, for instance a particle moving inside an harmonic plus a quartic term potential V(x) = — \x 2 + jX 4 . 
In this case the potential has two wells located at x = If two systems described by the variables x, y start to 

evolve within the same well (i.e. x, y > or x, y < 0) they will always tend to finish in the same final configuration 
because the synchronizing effect of the noise is, at very long times, the dominant effect. If they start in different wells 
the conclusion is also the same because there is always a finite probability that a strong fluctuation in the stochastic 
noise drives both particles in the same well. This conclusion, which holds for maximal cross-correlations JC = 1, can 
be generalized for any potential with a finite number of wells separated by finite energy barriers. There is always a 
finite probability that a fluctuation of the noise can take both particles into the same well and hereafter the distance 
D(t) between both configurations would tend exponentially fast to zero. Obviously this argument applies only for 
finite barriers, finite wells as well as finite temperature. At zero temperature the synchronizing effect of the stochastic 
noise is absent and the asymptotic distance may not vanish and show a non-trivial dependence on the initial distance. 

From the discussion above it follows that the role of the thermodynamic limit N — ► oo in the DS is also crucial. 
In this limit the height of the barriers or the number of wells (i.e. meta-stable states) may diverge. The first case 
happens, for instance, in the ferromagnetic Ising model where the time-reversal symmetry of the Hamiltonian is broken 
below T c . The second case is realized in spin- glass models where the number of meta-stable states is exponentially 
large with N . In both cases, a finite fluctuation of the stochastic noise (even with JC = I) may not synchronize the 
system and Doo can be a non trivial function of both the temperature and the initial distance. To be more precise 
Doq is defined as follows, 

= lim lim D(t) (7) 

Note that for N finite we expect limt_>oo D(t) — so it is crucial that the thermodynamic limit N — > oo is taken 
before the infinite-time limit. Taking the limits in the reverse order will result in that will always vanish at finite 
temperature. Note that this discussion applies only when the D = stationary solution is stable. This is indeed 
satisfied for JC = 1 but may be also fulfilled in more general situations with JC < 1 and D = still being an stationary 
stable solution. 

There is another property of the case JC — I which makes it particularly interesting. Up to now our discussion 
was limited to different mathematical properties of damage spreading transitions. But, what about their physical 
significance? Suppose we take two generic statistical systems described by the set of variables {xi,yi},i = 1, ..,N 
which evolve under the Langevin dynamics, 

^l=F i ({x(t)}) + r 1t ; ^ = F i ({y(t)})+r,' t (8) 



where, as before, ry, rf are white noises with cross-correlation JC which we will suppose is a generic function of the 



Hamming distance D. Let us suppose that the force derives from a potential Fi({x}) — — dv ^ x ^ . if we define the 



new variables z% — Xi — yi we may obtain, subtracting both equations in (g), 

= F t ({y(t) + z(t)}) Fi({y(t)}) + u t (t) (9) 

where v is a stochastic white noise of zero mean and variance 4T(I — JC(x,y)). The solution Zj = is a stationary 
solution of (|9|). A linear stability analysis around that solution yields the equation, 

-£=H i3 {{y{t)})z j {t) + v i {t) (10) 

where Hij (y) — — is the Hessian matrix evaluated at the point (yx, ■■, 2/jv)- The solution Zj = 0, correspond- 

ing to vanishing Hamming distance D = 0, is stable if the Hessian is negative definite. The presence of the stochastic 
noise v in ( |l0| ) decreases the stability of the D — solution. Because JC < I we conclude that the region where the 
D = stationary solution is maximally stable corresponds to the case when JC = 1 because v vanishes. So JC = 1 is 
the cross-correlation for which the stationary solution D = is maximally stable. For JC — I equation (|l^) is quite 
appealing and evidenciates the physical origin of the DS transition. An instability of the D = solution may appear 
when an eigenvalue of the Hessian matrix Hij vanishes. This corresponds to a saddle point of the potential land- 
scape V({y}). Due to ergodicity, the x,y systems sample all the possible configurations. So the asymptotic distance 
Doe is a direct measure of the stability of the D = solution along all possible configurations (weighted with their 
corresponding statistical Boltzmann weight). In other words, a instability in the solution = and the existence 
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of a DS transition is indication of the presence of saddle points in the potential landscape of systems x, y. If K, < 1 
the temperature of the DS transition will depend on the particular form of K, (actually we will see later, in the study 
of mode-coupling equations, that it depends through the value K'(D — 0)). Furthermore, due to the destabilizing 
effect of the noise the damage spreading transition will increase when K, decreases so JC — 1 yields the lowest damage 
spreading transition temperature among all possible cross-correlations IC(D) for which the D = solution is stable. 
In the presence of fixed points for the asymptotic distance other than D = the stationary solution = may 
become unstable because the noise v is too strong (similarly to what happens in the harmonic oscillator example) and 
damage spreading does not evidenciate anymore the existence of saddle points. In other words, saddle points may be 
observed only by studying the T dependence of the basin of attraction of the Q = stationary solution supposing it 
stable at very highptemperatures (free case). Remark that a similar argument has already been presented by Loreto, 
Serva and VulpianiE.il for systems described by a single variable x(t) in a potential field V(x). 

We have argued above that for Langevin dynamics the maximal cross-correlator K. = 1 is a special case, resulting 
in a simplification of the problem. We emphasize that in other dynamical systems it is unclear whether or not the 
maximal cross-correlator K — 1 plays the same role in the context of damage spreading. This is due to the complexity 
of cross-correlations. 

In the next section, we will analyze in detail the DS dynamics in the mode-coupling theory of glasses. As has been 
previously said these equations describe the relaxation processes and dynamics in glasses (in the under-cooled liquid 
regime) in the mean-field approximation and represent the dynamical behavior of systems with an exponentially large 
number of mcta-stable states. 



III. DAMAGE SPREADING IN MODE-COUPLING THEORY 



Mode-coupling theory describes relaxational processes in glasses. In the nutshell, mode-coupling theory corresponds 
to an exact re-summation of an infinite series of diagrams in the hydrodynamic theories. The kind of diagrams which 
the mode-coupling approximation selects are those which precisely survive in the mean-field limit of some realistic 
models. So, a way. to get mode-coupling equations is by considering exact dynamical theories for mean-field disordered 
spin-glass modelso. Spherical spins allow for an exact closure of the djaiamical equations in terms of correlation and 
response functions (as was shown by Crisanti, Horner and SommerscJ in the p-interactions spherical spin glass), 
leading to many analytical results. Although spherical spins are unrealistic (compared to Ising spins) they capture 
the essentials of the dynamics which is universally found in a large variety of models. Whereas for p = 2 the physical 
description of the model is quite simpldI9 the behavior turns_out to be much more interesting for p > 2 where an 
exponentially large number of meta-stable states are presenter. In this type of models meta-stability plays a very 
important role so, according to the arguments of the previous section, we expect to get interesting results for the DS 
transition. Forthcoming subsections analyze this transition in detail. A good review on the main results obtained in 
this model have been collected and reported by BarratE3. 



A. Random initial configuration 



This section is devoted to the study of the DS problem in the mode-coupling equations. This section describes 
some preliminary work already presented ino but here we present a more extended research of the problem to include 
asymmetry, different class of initial conditions as well as general cross-correlations of the noises. The simplest solvable 
model whose dynamics is described by the off-,equilibrium mode-coupling equations is the spherical p-spin glass model 
introduced by Crisanti, Horner and SommersE3. In this case, the configurations are described by N continuous spin 
variables {erf, 1 < i < N} which satisfy the spherical global constraint Yli=i °f = N. 

The Langevin dynamics of the model is given by, 

^i= Fl ({a})-^ l+Vl (11) 
where Fi is the force acting on the spin Oi due to the interaction with the rest of the spins, 

^ = 7^1)7 £ Jf^^-^ ■ (12) 

(i 2 ,i 3 ,...,i p ) 

The term fi in eq.(|l^) is a Lagrange multiplier which ensures that the spherical constraint is satisfied at all times 
and the noise 77 satisfies the fluctuation-dissipation relation (r]i(t)r)j(s)) — 2TS(t — s)5ij where (...) denotes the noise 
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average. The J*' 3 '"' p are quenched random variables with zero mean and variance p\/(2N p ^ 1 ). The interactions 
jvi, 3,-, p are S y mme f ;r i c un der the interchange of the super-indices «2j*3j ■•>%) but in the most general case may not 
be symmetric under the exchange of the subindex i with a generic super-index. So, for instance J i 2,l3, ^^ , p ^ j^'-^p _ 

Most of the studies undertaken in this model concentrate on the symmetric case where J* 2,l3 '"' 4p is symmetric under 
the permutation of all possible indices. This case is particularly interesting because there exists an energy function 
such that the force Fi derives from a Hamiltonian or potential function = — J2£ so there exists an stationary state 
described by a Boltzmann-Gibbs distribution. Due to the mean-field character of the model, the dynamical equations 
depend on the statistical properties of the force only through its correlations. OaJhe other hand, the statistical 
properties of the force Fi depend on the correlations of the J's. The simplest caseO corresponds to correlations of 
the type, 



t»2!*3j">»j> T»l>")t*-l)*lle+li"|i|> _ -P- f1Qt 

for every k. So if a = 1 we recover the symmetric case-while for a = we obtain the asymmetric case. Eq.(|l3|) 
implies the following statistical properties for the force FiEB, 

WMWjWF) = Sijf'iQ) + (1 ~ (14) 
where f(q) = q p /2. In the asymmetric case a = the forces are completely uncorrelated at different sites. Hence 



the set of equations (11) become uncorrelated and the problem can be partially solved. This particular case will be 
analyzed later. For a < 1 there does not exist an energy function TL which drives the system to thermal equilibrium 
and the fluctuation-dissipation theorem is not fulfilled. 

We define the overlap between two configurations of the spins cr, r by the relation Q = - X)i=i a i T i so * ne Hamming 
distance between these two configurations is, 

B = ^ (15) 

in such a way that identical configurations have zero distance and opposite configurations have maximal distance 
D = 1. Then we consider two copies of the system {ui,Ti} which evolve under the same statistical noise ( |Tl| ) with 
cross-correlation K, but with different initial conditions. We assume the cross-correlator to be a function of the 
Hamming distance D or the overlap Q. As a consequence, other choices for the cross-correlator (in general this could 
depend on both configurations C,C) may change our results obtained below for K,(Q) ^ 1. The major part, however, 
is concerned with K, — 1 (we will explain why) and will not suffer from this. In this section we restrict our attention 
to random initial configurations (i.e equilibrium configurations at infinite temperature) with initial overlap Q(0). The 
case of initial equilibrium configurations will be analyzed in the next subsection. The different set of correlation 
functions which describe the dynamics of the system are given by, 

N N 

C(t,s) = (l/iV)5>^K( s )) = (l/N^initMs)) (16) 

i=l i=l 

R(t,s) = (l/N)jr^ = (l/N)jr^- (17) 

i— 1 1 i—l % 

N 



Q(t,s) = (l/N)J2(<ri(t)n(s)) (18) 



i=l 



where < .. > denotes the average over dynamical histories and /if, hj are fields coupled to the spins Ui, respectively. 
In what follows we take the convention t > s. The previous correlation functions satisfy the boundary conditions, 
C(t,t) = l,R(s,t) — 0, lim t _ > ( s )+ R(t, s) — 1 while the two replica overlap Q(t,s) defines the equal time overlap 
Qd{t) — Q(t, t) which yields thaJJamming distance at equal times or damage D{t) through the relation (|l5|). Following 
standard functional methodsBt2l it is possible to write a closed set of equations for the previous correlation functions. 
Some details of the computation are shown in Appendix A. The final result is: 

<(i{t)C{t,s) = ^ [ duR(s, u)C p - 1 (t, u) + a P( - P ~ X) / duR(t, u)C(s, u)C p ~ 2 (t, u) (19) 



dt ^ w v ' ' 2J 
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dR ^ S) + fi(t)R(t, s) = S(t -s) + a P( - P ^ / tt)i2(«, s)C p - 2 (t, u) (20) 

at 2 Js 

0Q(f, s) + pffiQfo fi) = P I" dui?(S; u )gp-i( t) u ) + a pi - p - X ) / duj^ u )g( Mj s )C p ~ 2 (i, u) (21) 



at ' 2 

while the Lagrange multiplier and the diagonal correlation function Qdit) obey the equations, 



fi(t) =T+ PH±^k 11) f duR(t,u)C p - 1 (t,u) (22) 
2 Jo 

+ fi(t)Q d (t) = TK{Q d {t)) + | jf dufl(f, u)Q p - 1 (t, u) + a p(p ~ 1} jf dufl(t, u)Q(t, U )Cf- 2 (t, u) . (23) 

Note that the cross-correlation fC{Qd) only enters explicitly through the equation ([23] ) so it does not affect the 
evolution of a single replica. The whole set of equations is quite involved. For the correlation C and response 
functions R eqs.([19 20 22j) several results are known, in particular their behavior in the equilibrium regime (where 
time-translational invariance is satisfied and the fluctuation-dissipation theorem is obeyed) as well as in the non- 
stationary aging regimeEi 

In what follows we analyze different dynamical fixed points of eq.(|2^) and show the existence of a dynamical 
instability in the DS equations. 



B. Fixed-point analysis for a generic cross-correlation K. 



Different type of dynamical regimes may be distinguished depending on the cross-correlator JC. Our analysis is 
similar to that performed in section II for the simple harmonic oscillator. Different fixed points for the dynamical 
equations can be analyzed from equation (E3T). If the temperature T is very large then eq.(E3l) becomes, 



\^P-=T(K(Q d (t))-Qd(t)) (24) 



where we have used fx = T using cq.(|22]). Equation (24) can be exactly solved. The stationary solutions are given by 
)C(Q) = Q- In figure (Q) we analyze the different solutions for a generic JC. We find that there are different stationary 
solutions corresponding to all possible intersections of the two curves (Q and JC(Q)). A linear stability analysis of (g4) 
reveals that stationary solutions Q* are stable if tC'(Q*) < 1. A dynamical flow diagram can be constructed where 
the region of stability is indicated by different arrows. Regions where JC{Q) > Q satisfy > and regions where 
tC(Q) < Q satisfy < 0. So in this case one may depict a diagram of all possible high-temperature dynamical 
phases which separate regions with different fixed-point attractors. Damage spreading transitions will strongly depend 
on the type of cross-correlator. The case JC — 1 is shown in figure (Q) where there is a unique attractor at Q — 1 at 
very high-temperatures. This analysis of the different dynamical phases is valid only at very high temperatures. As 
soon as the temperature is finite and starts to decrease some of the stable fixed points may become unstable and other 
unstable points may become stable. The damage spreading transition correspond to the appearance of an instability 
in one of these high-temperature fixed points. As we will see below, the damage spreading transition temperature may 
be different for different fixed-points since it depends on the value of JC'(Q*) which may vary for different fixed-points 

or- 
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FIG. 1. Flow diagrams and fixed points for a generic correlator /C at infinite temperature. Q* and I are stable fixed points 
and and the fixed point between Q* and I are unstable. 

In what follows most of our discussion will concentrate on the particularly interesting case fC = 1 which has a 
unique fixed point at Q* = 1. Although the analysis may be extended to other fixed-points this case is also the 
most interesting according to our preceding discussion in section II. As we will check below, this case also defines the 
lowest damage spreading temperature To among all the possible cross-correlators K, for which the fixed-point Q* = 1 
is stable. 

A 




Q Q*=l 
FIG. 2. Flow diagram for K, = 1. 



1. Existence of Tq : lower and upper bound for a — 1 

A first glance to equations (|2l|) , (p3|) reveals that the overlap Q(t 7 s) and its diagonal part Qd{t) are coupled each 
other through the correlation C(t, s) and response function R(t, s). The trivial solution Q(t, s) — C(t, s) and Qd(t) = 1 
corresponds to the case where the initial conditions are the same, Qd(0) = 1 and the distance D(t) = for all times. 
This high-temperature fixed-point (hereafter we will denote it by HT) corresponds to = and is asymptotically 
reached by the dynamics for high enough temperatures under certain conditions of the cross-correlator K. (see the 
preceding discussion). In what follows we concentrate our attention in the case a = 1 where there is a stationary 
solution for C, R and /x corresponding to the equilibrium measure. Numerical integration of those equations (see later) 
reveals that the typical time needed to reach that solution grows if temperature decreases. At a given temperature 
(which we identify with To) there is an instability in the dynamical equations (pT|) , (|23|) and the asymptotic solution 
differs from the HT one. We did not succeed in finding an explicit expression for To but we have been able to show 
its existence and find a lower and upper bound for its value (see Appendix B): 

^2(1 -£'(1)) ~ T ° ~ V 2(1 -K'{1)) (25) 

Note that for the particular case p = 2 and /C'(l) = the inequality (^5|) yields Tq < 1. Taking into account that 
( p5| ) was derived under the assumption To > T s = 1 (i.e. we.-supposed we were in the high-temperature regime) this 
is not inconsistent with the result To = 1 derived by StariolcEl 

Note that both the lower and upper bound for To diverge when ZC'(f) = 1. This limiting value sets a condition 
on the possible cross-correlations JC{Q) where the Q = 1 solution is stable. Only those functions IC(Q) for which 
ZC'(l) < 1 are those for which Q = 1 is linearly stable at very high-temperatures. According to our discussion in 
section II the appearance of DS in this case is related to the presence of a divergent number of saddle points which 
mark the onset of a dynamical instability. Because 1C'(1) < then (according to (p5|)) the limiting case /C'(l) = (for 
instance, if IC(Q) — 1) sets the lowest value of the temperature To where DS appears. This is important because it 
means, that whatever correlator JC we consider (such that the solution Q = 1 is stable for high-enough temperatures) 

damage spreads below To = \J~^^-- Remark that in the general case p > 3 the dynamical instability temperature To 
stays well above any relevant critical temperature (T s or Td)- 

In the next section we discuss the behavior of the asymptotic distance as a function of temperature. For simplicity, 
our analysis is restricted to the case JC = 1 for which most of the numerical work has been done. We will see that 







-Dqo, for a given specification of the correlator /C, seems indeed to play the role of a dynamical order parameter in DS 
transitions. 



2. Numerical analysis 

In general it is too complicated to obtain an analytical solution of the set of equations (|l9|)~(p3|). We shall devote 
this section to a numerical study of the equations ([L9|)-(|2^) for the DS problem. Although in some particular cases 
an exact solution can be found (see below) this is not the general situation. 

First, one could investigate the long-time limit of Qd via a numerical integration of the set of equations jl9|)- 
(p3|). However, the CPU time and the memory needed to do this grows very fast with time because of the integrals 
occurring in the equations. Thus the spreading of damage at large times can only be obtained from the dynamical 
equations doing some extrapolations. This enlarges the error in the estimate of Z?oo, especially in cases where Qd(t) 
is a non-monotonic function of time. In figure (0) we show how the overlap Qd(t) depends on the initial condition. 
Although direct extrapolations from numerical data of the value of the asymptotic damage are difficult, the figure is 
not incompatible with an independence of on the initial condition. Another more powerful technique is necessary 
to corroborate this result. 



FIG. 3. Qd(t) for p = 3 (a = 1, /C = 1) at temperatures T = 0.1, 0.5 (from bottom to top at large times) for three different 
values of the initial overlap Qd(0) — —1, 0, 0.5 as a function of time. The continuous lines are the numerical integrations with 
time step At = 0.01. 

An alternative approach to obtain the long-time behavior of time dependent variables with high accuracy was 
introduced by Franz, Marinari and ParisiEj to study the long-time behavior of the energy. Here we extend their 
method to analyze the asymptotic behavior of D{t). In their method they first decompose the time dependent 
variables in a series expansion before extrapolating for large times with the help of Pade approximants. For the DS 
problem, it leads to a Taylor expansion of the correlation function C, the response function R and the overlap Q: 



oo oo 



oo oo 



oo oo 



(26) 



k=0 1=0 



k=0 1=0 



k=0 1=0 



where coo = foo = 1 since C(t, t) = 1 and lim t ^( s )+ R(t, s) — lEa. Moreover, fi(t) and the diagonal correlation function 
Qd{t) can be written as: 



fc=0 



oo k 



i(t) = j2^t k , Q d (t) = EE 9 ( fc -o'* fe ' 



(27) 



k=0 1=0 



with fiQ — T andgoo is a parameter identical to the value of Qd at t = 0. Assuming always t > s, the dynamical 
equations (|l9|)~(p3|) can be transformed into recurrence relations for the coefficients of the expansion. To this end 
one first substitutes (^) and ( p7j ) into ([l9|)-(p3|) under the constraint t > s, then calculate the integrals and finally 
rearrange terms. 
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Numerically, the coefficients of the expansion are now readily obtained. In case p = 3, the first 80 coefficients of the 
expansions can be computed on a RISC workstation in a few hours. However, for larger values of p the computational 
effort is larger. 

To ensure high accuracy of the asymptotic extrapolation, one needs a large radius of convergence of the series 
expansion. A good method to enlarge the radius of convergence of a series expansion is using Pade approximants. In 
this method one introduces two polynomials U m (t) and Vk(t) of degree at most m and k, respectively. The goal is 
to choose U m (t) and Vk(t) for given m and k such that Qd{t) and U m (t)/Vk(t) are equal at t = and have as many 
equal derivatives as possible at t = 0. 

The computations have been performed for the symmetric case a = 1 and JC = 1. Moreover, three different initial 
conditions have been considered: a) anti-correlated random initial conditions with Qd{0) = — 1, b) uncorrelated 
random initial conditions with Qd{0) — and c) partially correlated random initial conditions with Qd{0) = 0.5. 

To check that the extrapolations using the Pade approximants are correct, the Pade series have been compared 
with numerical integrations of the dynamical equations. Indeed, the Pade series and the numerical integration fit 
closelyM 

The estimate for is obtained by division of the highest order coefficients of P m (t) and Qfe(i), i.e., by a m /bk- 
Moreover, an asymptotic estimate can be obtained assuming a power law decay of the equal time overlap: Qd(t) = 
Qd(oo) + At 1 . The analysis of £>oo suffers in some cases from a small radius of convergence (even with Pade) as well 
as from the presence of poles in the Pade expansion. The results are displayed for p — 3 in figure ^ and for p — 4 in 
figure ^ for cases a), b) and c) as a function of the temperature. Let us remark that a lower number of coefficients 
in the Taylor expansion in the case p = 4 with respect to p = 3 leads to a less accurate estimate of the asymptotic 
distance. 
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FIG. 5. Asymptotic distance D^ for p = 4 (a = 1,JC = 1) obtained from the Pade analysis of the series expansions for 
different initial conditions Do = 1 (crosses), Do = 0.5 (triangles), Do = 0.25 (stars). Typical error bars are shown for the case 
Do = 0.5. 

Inspection of figures |] and |E| reveals that the dynamical transition To is in the predicted regime eq. ( |25| ) . It can be 
estimated more accurately from the relaxation time T re i ax associated to the decay of the distance D(t) to zero. Starting 
from high temperatures, we assume the relaxation time to diverge at To according to a power law: r ro i ax ~ (T— Xb)~ 7 . 
We thus have found T (p = 3) = 1.04 ± 0.02 with 7 ~ 1.1 ± 0.1 and T (p = 4) = 1.13 ± 0.02 with 7 ~ 1.1 ± 0.1. 

We conclude that for all temperatures, both p = 3 and p = 4, the asymptotic distance is independent of the initial 
distance. This is ip, contrast with the case p = 2, where dependence of the initial distance is found for the low 
temperature regionta We must say that we have obtained the same results, as in cases p = 3, 4, for a model which 
is a combination of the p = 2 and p = 4 spherical spin-glass modeo. For a certain range of parameters, this model 
is known to have a continuous phase transition with continuous replica symmetry breaking and without collapse of 
the configurational entropy. So the first-order character of the spin-glass transition found in our model for p > 3 is 
not essential for the appearance of the DS transition. Still, that modec3, is also characterized by the presence of an 
exponentially large number of meta-stable states. From the point of view of the form of the dynamical equations, 
the fact that To is present for p > 2 as well as in a model which combines p = 2 and p — 4 is a consequence of the 
non-linearity in the coupling between the damage Qd(t) and the two-time correlation function Q(t, s) which occurs 
for all p > 2. From the physical point of view, this independence of on the initial distance is quite appealing. 
In general, one would expect that p = 4 is quite similar to p = 2 due to the presence of the time-reversal symmetry. 
The fact that the damage does not feel this symmetry for p > 2 means that the separation of dynamic trajectories 
does not occur in the borders or maxima which separate equilibrium states but within saddle points of the phase 
space which divide configurations separated by finite energy barriers. This is supported by the fact that the transition 
occurs at a temperature much higher than T s and, as we will see in the following section, from the fact that it happens 
starting already from an equilibrium configuration. The asymptotic value Z?oo can, on the basis of our computations, 
be regarded as an order parameter for the transition at Tq. Although Doo and the transition Tq itself do depend 
on the specific choice of the correlator JC it is interesting that the asymptotic state does not depend on the initial 
distance. For a better understanding of the physical origin of this transition we shall consider the case of equilibrium 
initial conditions in the next section. 



C. Equilibrium initial condition 

We have seen in the previous section that the asymptotic distance Doo is a non-trivial function of the temperature 
which is finite below Tq and vanishes above T). The relevance of the existence of the meta-stable states has been 
already pointed out in previous sections specially when the cross-correlator is maximal JC = 1. The fact that the DS 
transition exists suggests that the nature of this phase transition is related to the corrugated properties of the free 
energy landscape. To check this result it is convenient to investigate the DS transition starting from an equilibrium 
configuration. In this case the system starts from a stationary state and remains there forever. At high temperatures 
this state is the paramagnetic so in this case the DS is a direct check pi the ruggedness of the paramagnetic state. In 
fact, we will find that if we start from an initial equilibrium conditioned, then the DS transition persists and actually 
coincides with the previous Tq found for a random initial configuration. This reinforces the idea of Doo as a dynamical 
order parameter for the DS transition for a given choice of the correlator JC. Again, for sake of simplicity, we restrict 
our analysis here to the case tC = 1. 

The analysis of the dynamical equations for an equilibrium initial condition follows the same steps as for the random 
case but now we must impose a Gibbs distribution for the configurations a and r at time 0. Nevertheless, there is 
a point which must now be considered. Let us take T > T s (i.e we will suppose equilibrium configurations in the 
paramagnetic phase). Suppose that p is odd and we take an equilibrium configuration at temperature T. To impose 
Qo = — 1 or Do = 1 we must take <Ji = — Tj for all i. Because the energy is an odd function of the spin variables we 
have E({a}) = —E({t}). If the equilibrium energy is not zero (this happens everywhere except at (3 = 0) we cannot 
put both configurations at equilibrium at the same temperature (because the temperature uniquely determines the 
value of the equilibrium energy). Then, if both initial conditions a and r are equilibrium initial configurations their 
overlap Q(0) must vanish. Actually, for T > T s two equilibrium configurations do have overlap zero with probability 
1 and overlap different from zero with probability exponentially small with N. So if we take the thermodynamic limit 
before the infinite time limit it is clear that we must start with zero initial overlap. To be more precise, the probability 
that two equilibrium configurations {cr}, {r} in the paramagnetic phase do have overlap q = 53i=i <T i T i is given by, 

P(q) ~ exp(-NPf(q)) (28) 
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where f(q) is the free energy cost to find a correlation q between the configurations. Clearly, because q = 
corresponds to the equilibrium value in the paramagnetic phase, f(q) has a minimum at q = so we can write 
f(q) = constant + q 2 /(2(3xsa), 



, Nq 2 . 
P{q) ~ exp(~- ) 

*XSG 



(29) 



with xsg — N < q 2 > is the spin-glass susceptibility. Above T s the xsg 1S finite and the probability to have q =^ 
is exponentially small with N. Now the cost in free energy f(q) has two parts, a cost in energy u(q) and a cost in 
entropy s(q) — (3(u(q) — /(?))■ The cost in energy vanishes at infinite temperature and the full cost of f(q) is due 
to the entropy. So only at infinite temperature (i.e random initial configurations, the case considered in preceding 
sections) we can impose an initial condition with initial non-zero overlap. 

The equivaleat of the dynamical equations eqs.([H) 20 ^) can be easily obtained for T > T s in the replica symmetric 
approximation^. The correlation function C(t,s) and the response function R(t,s) are time-translational invariant 
and satisfy the fluctuation-dissipation theorem TR(t) = 8(f). The C(t) satisfies the equation, 



duCP -\t-u)^±=0 



(30) 



with C(0) = 1. The two-times overlap satisfies the equation, 

dQ(t,s) , (rr ^PP, n( . \ PP f s . nP -i u , dC( S ~u) p0 f* dCP-Ht-u) 

- {T + —)Q(t,s) - — I duQ p (t,u) — — / duQ{s,u) = 



dt 



du 



(31) 



with the initial condition Q(t, 0) = Q(0, t) — 0. The diagonal part, Qd(t) — Q(t, t) is given by, 



-T 



Pi 
2 



Oil 



pP 

2 



, ,dCP- l {t-u) n 
duQ{t,u) ^ '- = 



(32) 



T 



We have looked for a time- 



with the initial condition Qd{0) — 0. Now we are in equilibrium so fi(t) 

translational invariant solution for Q(t, s) (i.e a solution of the type Q(t, s) = Qd(s)Q{t — s) for t > s) but we have 
not found it (even for p = 2). Our numerical results suggest that such an asymptotic solution does not exist. 

Using as before a series expansion in the time-dependent variables and Pade approximants, we have estimated the 
asymptotic distance for equilibrium initial conditions, i.e., for Q(0) = 0. The results are displayed in figure ^ for 
different temperatures. The divergence of the relaxation time leads to To(p = 3) = 1.01 ± 0.04 with 7 = 1.4 ± 0.3 
which indicates that To coincides with the result obtained for random initial conditions. This supports the idea that 
the transition at Tq is of a dynamical nature and unrelated to the existence of a thermodynamic phase transition. 
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FIG. 6. Asymptotic distance Doo as a function of the temperature, starting from an equilibrium configuration, for p = 3 
(a = 1,K= 1). 
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D. The non-symmetric q/1 case 



As we saw in the previous section, one of the most interesting results concerning the DS transition is the fact 
that D(po) is a non-trivial quantity which does not depend on the initial condition (and depends only on the cross 
correlation /C). In equilibrium thermodynamics, this is one of the features of order parameters which separate different 
equilibrium phases. In the stationary state, when fluctuation-dissipation theorem is obeyed, the order parameter is a 
a quantity which characterizes the equilibrium state and (in the absence of ergodicity breaking) does not depend on 
the initial condition. 

In order to present a convincing proof on this result we have investigated the general non-symmetric case a 1. A 
difficulty inherent to the extrapolations made from figures 1 and 2 is the fact that, below To, the convergence of the 
distance D(t) towards its asymptotic value is very slow (a power law in time). Consequently, both numerically or 
using the Pade method it is very difficult to extrapolate to the asymptotic value. As the asymmetry of the interactions 
eq.([f|) is turned on (i.e if a < 1) the relaxation of the system to the stationary state turns out to be faster. Actually 
for the asymmetric case (a = 0) or the antisymmetric case (a = — jzj, see later) the relaxation of the distance D(t) 
is nearly exponential. We have no reasons to suppose that the independence of the asymptotic value on the initial 
distance is a dependent. Our analysis for a < 1 suggests that the independence of D(t) on the initial distance Dq 
holds for all generic values of a. 

It is important to note that, for a < 1, there is no equilibrium stationary state and the fluctuation-dissipation 
theorem is not satisfied. Still we expect, for p > 2, the DS transition to survive for any quantity of finite asymmetry. 
The reason is that the DS transition (even for a = 1) is inherently a non-equilibrium transition so the effect of the 
asymmetry may not change the character of that transition. i_ . 

The non-symmetric case for p — 2 has been already considered by Crisanti and Sompolinsky ten years agcc3. By 
assuming that, in the stationary state, the correlation and response functions are time-translational invariant they 
succeeded in showing that the thermodynamic transition T s = 1 for a = 1 turned out to be unstable against the 
asymmetry for any value of a < 1. They also derived the result fi(oo) = \/l + T 2 for a — 0. 

In the next paragraphs we study the cases a — 0, starting from a random initial configuration in some detail. 

Unless stated we will consider the case JC = 1. 



1. The asymmetric case a = 

The case a = is quite interesting. The equation for the response function (|2^) simplifies considerably, 

8R(t,s) 



function becomes, 



dt 



+ fj,(t)R(t,s) = 6{t- s) (33) 



which can be readily written, using R(t,s) = ^fey with z(t) = exp(fy yn(t')di') . The equation for the correlation 



dC(t, s) 



(34) 



dt ^ v ' y ' ' 2J z(s) 
Define the new function A(t, s) = z(t)C(t, s)z(s). In terms of this new function eq.(p4|) is, 



dt 2 J Z P- 2 (u)zP- 2 (t) 



From eq.(22) it is easy to derive a equation for z(t), 



2 dt 2 J Z P- 2 {u)zP- 2 (t) 

Equations eq.(p5|),(|3^) form a closed set of equations which can be solved with the initial conditions Aft, t) — z 2 (t), 
z(0) = 1. Once this set of equations is solved one can also find a solution for the overlap Q(t, s) in cq.(|2l|). Again, we 
define B(t, s) = z(t)Q(t, s)z(s) which satisfies the equation, 

dB(t, S ) P r du BP-\t, U ) (37) 



dt 2 J Z P- 2 (u)zP- 2 (t) 
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and the equal-times overlap eq.(|23J) b(t) 

1 db(t) 

2 dt 



B(t, t) satisfies the equation, 

T Z \t) + p -fdu ^Mn 



(38) 



with b(0) = Qd{0)- Note that this set of equations is quite involved for p > 2. Only for p = 2 they dramatically 
simplify (the case considered by Crisanti and Sompolinsky) and become linear. For general p the previous equations 
are non-linear. We have not succeeded in finding the asymptotic solution of these equations although we have guessed 
the results from the numerical results. We find that the DS transition is still present at finite temperature for p > 2. 
The analytical expression for T c is given by 



(39) 



Note that T c coincides precisely with the lower bound previously derived in eq.(p5|). The asymptotic value of fj,(t) 
is given by ju(oo) = \/l + T 2 and is p- independent. The asymptotic distance for p = 3, 4 is given byEZl, 



£>oo = 



1 



T-(l-T c )m S 



(40) 



A full theoretical derivation of this results remains an interesting open problem. In what follows we compare our 
results obtained by numerical integrations with time step At = 0.01 with these theoretical guesses. In figure we 
show the overlap Qd(t) as a function of time for different temperatures below T c for p — 3. Note that the asymptotic 
value clearly does not depend on the initial condition. The horizontal dotted lines correspond to the asymptotic value 
eq.(^). This figure unambiguously demonstrates that the asymptotic distance does not depend on the value of the 
initial overlap Qd(0). In figure || we show /i(t) for p = 3, 4 compared with the theoretical prediction /i(oo) = \/l + T 2 . 




FIG. 7. Qd{i) (p = 3, a = 0, K. — 1) with Qd(0) = 0,0.5, —1 as a function of time for temperatures T = 0.1,0.3,0.5 (from 
bottom to top at large times). The horizontal dotted lines are the theoretical guesses. 
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FIG. 8. (j,(t) for a = 0, /C = 1 (p = 3 solid, p = 4 dashed) with Qd(0) = — 1 as a function of time for temperatures 
T — 0.1, 0.3, 0.5, 0.7 (from bottom to top). The horizontal line is the theoretical prediction. Note that ^t(O) = T. 



2. The antisymmetric case a — — -ij 



This is an extremal case where J ix ' v J t ~ ' '"' p is maximally negative. Physically this means that the 
force felt by an spin i due to a multiplet M. of p — 1 spins is as much opposite as possible as the force which feels 
another spin contained in that multiplet due to the action of another multiplet M' of p — 1 spins constructed from 
the rest of the p — 2 spins in the previous multiplet M. plus the spin i. In the particular case p = 2 this can be easily 
achieved making j( — —jf which corresponds, according to eq. (|l3|), to a = — 1. But in general, a can never be 
equal to —1 for p > 2. Take p — 3 and three couplings Jj fc , . It is clear that if jf k = — J] k and J] k = — j£ J 
then Jj fc 7^ ~Jk ■ It can be easily shown that the minimum value for a is given by a — — — ^j. 

Interestingly, this case can be exactly solved for the correlation and response function. Although it turns out to 
be quite difficult to solve for the overlap function we will analyze here a general correlator JC. For a — — g~^n e 1 s - 
( |l9| - |2"3| ) considerably simplify because /z(t) in eq.(p2|) is time independent. Because the initial configuration was taken 
random at time this means that the stationary state follows completely random configurations. 

To solve the equations define the following correlators, 

C(t,s) = c(t,s)exp(-T(\t- s|)); R(t,s) = r(t,s)exp(-T(t- s)) t>s (41) 

In this case the dynamical equations for C and R simplify. A particular solution for c(t, s) and r(t, s) can be found 
which simplifies dramatically the dynamical equations. This solution is given by c(t,s) — r(t,s) = f(t — s). This 
time-translational invariant solution is consistent with all dynamical equations for all times. The final closed equation 
for f(t) is given by, 



dt 



= - V - exp {-Tip - 2)t) f dup-\t - u)f{u) exp {-T(p - 2)u) (42) 
1 Jo 



with the initial condition /(0) = 1. But the equation for the overlap Q{t,s) is more complicated and cannot be 
reduced to a time-translational invariant solution. Writing Q(t,s) — q(t, s) exp (— T\t — s\) (with q(t,t) = Qd(t) we 
obtain the following equations, 



dq(t,s) = p 
dt 2 



duexp(-2T(s - «))(/(« - w)<f ^(M) - f p ^(t - u)q(s,u)) (43) 
duf p ~ l {t - u)q(u, s) exp (-T(p - 2){t - u)) (44) 



' " ( '- !ih ' T{Q d {t)-K,{Q d )) = l / rfuexp(-2T(t- U ))(/(i- U )^- 1 (t, U )-/f- 1 (i- U ) 9 (i, U )) . (45) 



ii 



2 dt ' 2 

A time-translational invariant solution for q(t y s) does not exist for all times (contrarily to what happens for C and 



R) because the first integral in the right hand side of eq. (44) does not vanish 



Previous equations are solvable for p — 2 and K, = 1. In this case we get, 

df(t) 



;)i , duf(t-u)f(u) (46) 

Q d (t) = 1 - (1 - Q d (0)) exp (-2Tt) (47) 



and there is no DS transition for p — 2 as expected (i.e. To = 0). For a generic cross-correlation K. let us note 
that the equation (|4^) reduces to equation (|24|) so there will be different asymptotic values depending on the value 
of K.'(Q*) at the different set of fixed points Q* = JC(Q*). It is notorious that the case p = 2, a = — 1 at finite 
temperature reduces to the infinite temperature case for any a. This is closely related to the fact that the stationary 
solution in this case coincides with the random initial configuration although this is not true anymore for p > 2. For 
p > 2 the lowest Tq temperature for maximal cross-correlation K, = 1 becomes finite and this is due to the rugged 
structure of the force landscape. 

Equations for the damage for p > 3 are difficult to solve. We have not succeeded to obtain an analytical expression 
for the asymptotic values as well as for T c and . Numerical integrations of the equations reveal that the transition 
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persists at finite temperatures for p > 2. Figure [| shows the overlap Qd{t) as a function of time for different 
temperatures for p = 3, JC = 1. In this case the transition is located at To ~ 0.595 ± 0.005. Note that again relaxation 
to the stationary state is faster than in the case a = 1 and Z?oo is again independent of D(0). 
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FIG. 9. Q d (t) (p = 3,a= -1/2, K = 1) with Q d (0) = 0,0.5, -1 as a function of time for temperatures T = 0.1,0.3,0.5 (from 
bottom to top). 

These results are quite appealing. Here we find a DS transition in the presence of a time-translational invariant 
solution for C and R, i.e when the system starts already in the stationary state. This is in agreement with the results 
of previous subsection III .B for a = 1 where a DS transition was found (at the same temperature than starting from 
random initial conditions) when the system already started in the stationary state. Figure [l^ summarizes our results. 
We show the T — a phase diagram of the DS transition for p = 3. 
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FIG. 10. T — a phase diagram for p = 3, K = 1. The line is a guide to the eye. 

Let us remark our final conclusion for this section. A DS transition is present for all models with p > 2 above the 
TAP temperature where an exponentially large number of states appear. For a given choice of JC such that Q = 1 
is stable at infinite temperature, this transition has the following interesting properties: a) The asymptotic distance 
is independent of the initial distance; b) Is also independent of the type of initial configuration; c) Is stable in the 
presence of asymmetry, (but is unstable for p = 2 in agreement with results derived by Crisanti and Sompolinsky). 
This suggests that Doo has some of the crucial properties to be a good dynamical order parameter. The correlator 
Q(t, s) is not time-translational invariant in the time scale in which the correlation and the response are. 

As we said previously we expect the properties of this transition to strongly depend on the type of dynamics through 
the cross-correlator JC. In the next section, we will discuss discrete (Monte Carlo) dynamics, the case in which different 
algorithms correspond to different cross-correlators. 
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IV. DAMAGE SPREADING IN DISCRETE GLASSY MODELS 



Up to now we have considered the DS problem in case of a dynamics continuous in time such as Langevin dynamics. 
Here we want to investigate damage spreading and in particular the existence of To for discrete dynamics such as 
Monte Carlo algorithms. 

The analogue of the functions K. for continuous dynamics are the different algorithms used in Monte Carlo dynamics 
for discrete dynamics. So different Monte Carlo algorithms determine different types of correlations between the noises. 
As for the cross-correlator K, this implies that the algorithms determine the structure of the high-temperature fixed 
points. Their instabilities determine the subsequent low-temperature behavior. As we have already commented in 
the introduction we therefore do not expect the DS transition to be universal and the results of this section aim to 
be compared with the results already obtained for the Langevin dynamics in a continuous system. 

One of the essential ingredients for the DS transition is the presence of two competing effects: a synchronizing 
force (the stochastic noise i] in the Langevin dynamics) and a landscape-dependent force which pulls configurations 
apart one from each other into different directions. In the case of a discrete (Monte Carlo) dynamics the equivalent 
role of the stochastic noise is played by the set of random numbers generated during the Monte Carlo updates. Now, 
the random number in the Monte Carlo algorithm (uniformly chosen between and 1) determines the probability 
of a move depending also on the configuration of the system. This last dependence corresponds to the role played 
by the cross-correlator K, in the Langevin case where the two noises i], rf may be different depending on the value of 
IC(C, C). Roughly speaking the Metropolis algorithm for Monte Carlo dynamics corresponds to the case K(Q) = Q 
for Langevin. It is easy to check that, at infinite temperature, the fixed-points in both dynamics are the same. 

For continuous (Langevin) dynamics we had the freedom to choose the maximal cross-correlation JC — 1. For 
discrete dynamics, however, this is not the case. There are several well known algorithms in the Monte Carlo 
approach according to which updating rule they use. For instance, Metropolis, Glauber, or heat-bath. Among these, 
the last one is the only one which has a unique fixed point Q* = 1 at infinite temperature. So, heat-bath dynamics is 
the closest case (but different) to the JC = 1 of Langevin dynamics. Here, our numerical investigation will focus on this 
type of discrete dynamics. Let us note that the other algorithms may show different behavior (due to the presence of 
other infinite-temperature fixed points) and consequently also different DS transitions. This non-universalito-,of the 
DS transition (as in our previous analysis of the Langevin case) has received some attention in the literatureEa. 



A. Damage spreading in the Bernasconi model 

Here we will analyze the Monte Carlo dynamics with the heat-bath algorithm for the Bernasconi modelH. This is 
a long-range interaction model without disorder which is known to have a glassy behavior being in the universality 
class of spin-glass models with one step of replica symmetry breakingEi Consequently, its dynamical behavior is the 
same as predicted by the mode-coupling theory. _ . 

The Bernasconi model (for simplicity we will consider the closed model, seeCJ for more details) consists of N Ising 
spins <7j = ±1 in a one dimensional chain interacting through a long-range four-spin interaction. It is defined by the 
following Hamiltonian, 

1 N ~ 1 

H =nJ2 C " ( 48 ) 
k=l 

with Cfc = y~lj—i Vj< J j+k and where we take periodic boundary conditions <Xj = Ui+N- 

In this model there are particular values of TV for which the ground state is exactly knownll. The interest of this 
model is that it behaves like a disordered spin glass .in. the-absence of explicit quenched disorder in the Hamiltonian. 
Apparently, disorder is self-induced by the dynamicscHcj c3. This means that dynamics itself generates slow evolving 
variables which effectively act as quenched disordered fields. This model is characterized by three temperatures- 
a melting crystal-liquid first-order transformation temperature Tm, a dynamical transition temperature Td ~ 0.5C3 
where the relaxation time diverges and ergodicity breaks and, finally, a static (or glass) transition temperature 
T s ~ 0.25 where replica symmetry breaks and the configurational entropy collapses (this is the ideal glass transition 
predicted in the AGM theory). 

In the heat-bath algorithm, to go from a configuration {<Xj(t)}, at a given time t, to a configuration {<Ji(t + At)} 
at the next time-step t + At with At = 1/N, a spin tjk is chosen at random among the iV spins to be updated. The 
probability to put the spin up or down is decided according to the intensity of the local field acting on that spin. 
More precisely, if we write the Hamiltonian ( f48| ) in terms of the local field H — — hk &k then the probability of 
putting the spin up (pk = 1) or down (ok = —1) at time t + At is given by, 
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P(a k {t + At) = a) = - + - tanh (/3h k (t)a) 



(49) 



with a = ±1 and h k (t) is the local field acting on the spin k at time t. Note that the probability ( f49| ) does only 
depend on the local field acting on the spin k and not on the actual value of that spin at time t. Then, a random 
number z(t) with a uniform distribution between and 1 can be introduced and spins are sequentially updated 
according to the dynamical rule 



a k (t + At) = sign 



1 



tanh (/3h k (t)) - z{t) 



(50) 



With this rule ( p0| ) we have studied numerically the damage spreading of three different initial conditions, as in the 
p-spin model: a) anti-correlated random initial conditions with D(0) — 1, b) uncorrelated random initial conditions 
with D(0) = 0.5 and c) partially correlated random initial conditions with D(0) = 0.1. For each of these cases, the 
distance D{t) is computed up to 100000 and 10000 Monte-Carlo time-steps for N = 1000 and N = 5000 respectively. 
To analyze the data, the logarithmic time with base a = 1.1 is considered. Moreover, the data is averaged in intervals 
of the form (a k , a k+1 — 1) with k a positive integer. For T = 0.3 the evolution of D(t) is plotted in figure 
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FIG. 11. The distance D averaged per interval (a k ,a k+1 — 1) as a function of logarithmic time with base a = 1.1 for 
temperature T — 0.3. The upper two curves are the result of initial condition D(0) = 1 for N = 1000 (boxes) and iV = 5000 
(diamonds). In the same manner, the middle (lower) curves are the result of -D(O) = 0.5 (-D(O) = 0.1). 

To obtain the asymptotic value -Dqo from figures like [ll], a power law fit for low temperatures is used. For high 
temperatures an exponential fit up to 2000 Monte-Carlo steps in real time is used. The results are displayed in 
figure ^2|. 
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FIG. 12. The asymptotic distance as a function of temperature for three different initial conditions: D(0) = 1 with 
N = 5000(1000) as closed(open) triangles, D(Q) = 0.5 with N = 5000(1000) as closed(open) circles and D(Q) = 0.1 with 
TV = 5000(1000) as closed(open) diamonds. 
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One observes that the distance does not vanish for any finite temperature. Moreover, figure [L2| indicates the 
existence of a temperature T\ above which the asymptotic distance is independent of the initial distance. Below this 
temperature Ti, Z?oo does seem to be dependent of the initial distance. This dependence is supported by a numerical 
extrapolation which could well fail when going to enormously large time scales. Still, what we certainly find is the 
appearance of a .dynamical transition temperature T± ~ 0.5 in very good agreement with the transition Td where 



The behavior we find here, when compared to the previous Langevin analysis for the p-spin model, may appear 
quite different. But a careful analysis reveals that this is not the case. If we consider that /C, for the Langevin case, 
is a generic function which may depend on the overlap as well as on the temperature T = 1/(3, we may then imagine 
a situation such as depicted in figure ( [l3|) where the infinite-temperature fixed-point Q = 1 becomes unstable as soon 
as (3 is finite. In this case the asymptotic distance would be a non trivial function of (3 and the damage-spreading 
transition could well happen at the usual dynamical transition Td where ergodicity is broken. The dependence of 
the asymptotic overlap on the initial value could be a consequence of the presence of different fixed points at low 
temperatures. 

In the most general case one could imagine a scenario with three possible different regimes. A high-temperature 
regime T > Tq where = independently of the initial distance Dq. An intermediate regime T± < T < Tq where 
Deo — £'co(7 1 ) is not zero but independent of the initial distance (this regime would correspond to the appearance 
of temperature dependent fixed point for (3 finite as depicted in figure (ft3|)). And finally a low-temperature regime 
T < T\ where — D^iT, D ) depends on both temperatura-and initial distance. The results we find for the 
Bernasconi model are the same as found by Derrida and Weisbucto for the Shcrrington-Kirkpatrick model. The fact 
that in this model T — > oo is related to the infinite-range character of the interaction. Actually, DerridaEj has found 
numerical evidence that for finite-dimensional spin glasses there exists a range of temperatures where the asymptotic 
distance vanishes and To is finite. The dependence of the asymptotic overlap on the initial condition found here and 
inlij below T\ could well be an artifact of the large-time extrapolation where the time window simulated and the size 
as well are not sufficiently large. Unfortunately, it is not easy to simulate very large times and sizes in infinite-ranged 
models like the present one. 



The results of this section show that, the DS transition is very close (and probably coincides) with the dynamical 
transition temperature Td (below which the system never attains equilibrium and ergodicity is broken). Nevertheless, 
in this case the asymptotic damage below 7\ apparently depends on the initial condition (and probably on the type of 
initial condition as well) although such a firm conclusion needs more understanding of damage-spreading transitions 
for generic updating rules. 



In this paper we have studied the problem of damage spreading in the mode-coupling theory of glasses. Mode- 
coupling theory is well known to describe relaxation processes in glasses in the mean-field approximation. A simple 
way to obtain the mode-coupling equations is analytically solving the dynamics of multipsin interaction spherical spin- 
glass models. These models are characterized by the presence of a huge number of meta-stable states (exponential 
with the system size) which appear at a temperature Ttap higher than Td (where ergodicity breaks) and T s (where 
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FIG. 13. Possible flow diagram at different temperatures for heat-bath dynamics 



V. CONCLUSIONS 
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replica symmetry breaks) . Because the phase space in this class of models is characterized by a extremely rugged and 
complex free energy landscape they are good candidates to study the landscape properties using techniques taken 
from dynamical systems. 

A very interesting technique which is able to probe the topological features of the phase space landscape is damage 
spreading. This consists in the study of the distance between the configurations of two stochastic systems submitted 
to the same realization of the stochastic noise but differing in the initial conditions. By the same realization of the 
stochastic noise we mean noises that are statistically identical although generally correlated through a function fC(Q) 
which satisfies the condition JC(Q) < JC(Q = 1) = 1. In general, any choice for the correlator JC alters the results. For 
Langevin dynamics, we have shown that interesting results appear for the case JC = 1. In that case, both noises are 
identical for the two copies independently of their configurations. This yields the lowest damage spreading transition 
for which the Q* = 1 fixed-point is stable at high-temperatures. Depending on the value of JC'{Q = 1) one finds a 
different damage spreading transition temperature up to the limiting case JC'(Q = 1) = 1 (see eq.(p5|)) where To = oo 
and the fixed-point Q* = 1 becomes unstable. Whether this holds for other type of dynamics is not studied and 
remains unclear. 

An exhaustive study has been done for the case JC = 1 although similar results are obtained for any JC for which 
the solution Q = 1 is asymptotically stable. In this case, through functional methods and using the Pade series 
expansion method (to make safe extrapolations in the asymptotic long-time limit), we have shown the existence of a 
damage spreading transition To in general mode-coupling equations with any degree of asymmetry in the interactions. 
In particular, in case of symmetric interactions (where detailed balance holds) we have found evidence for a damage 
spreading transition at a finite temperature. This transition occurs at temperatures To higher than Ttap, this last one 
being the temperature where an exponentially large number of meta-stable states start to appear. The transition is 
characterized by a dynamical order parameter which is the asymptotic distance between the two evolving replicas. 
Interestingly, has the good properties of order parameters being able to distinguish different dynamical phases 
(in our case, there are two possible phases depending whether D^, vanishes or not). These properties are: a) is 
independent on the initial distance Dq for a given class of initial conditions, b) -Doo does not depend also on the class 
of initial condition (whether they are random or thermalized) c) The DS transition is stable against the inclusion of 
asymmetry in the interactions (i.e. against the violation of detailed balance in the dynamics). Furthermore, regarding 
the mode-coupling equations with asymmetry we have obtained some exact results for the asymmetric case a = 
and exactly solved the correlation and the response function for the antisymmetric case a = — which interestingly 
turns out to be time-translational invariant. The existence of DS transition in this case reveals that this transition 
already appears when the system is timc-translational invariant. 

We insist on the fact that the precise value of the damage spreading temperature To as well as the asymptotic distance 
Doo both depend on the correlator JC considered. This fact expresses the non-universal character of this transition 
where the cross-correlator JC plays the equivalent role of a stochastic noise for the dynamical order parameter Qd(t). 
Different functions JC imply different dynamical phase transitions so their physical significance must be appropriately 
interpreted. In this direction we have tried to interpret our results in terms of saddle points in phase space for those 
cases where Q* = 1 is an asymptotically stable solution at high temperatures. The essentials of the argument where 
given in section II where it was shown that for close enough initial conditions an instability in the Q = 1 high- 
temperature solution is due to the presence of vanishing modes in the Hessian of the potential function. This signals 
the presence of saddle points in the free energy landscape, of which an infinite amount yields a DS transition. This 
certainly happens at Ttap but we have found that the transition occurs already at a temperature To much above that 
temperature. How can we reconcile our results with that? At Ttap the number of meta-stable states is exponentially 
large with N. Probably for the DS transition to occur it is only necessary that this number be a big (for instance 
, a power of N) but not exponentially large. Consequently, a divergent number of meta-stable states and saddle 
points (but not exponentially large with N) should be enough to make the DS transition appear. Unfortunately, the 
analytical computation of the temperature To by counting the number of meta-stable states is not so direct because 
such a calculation involves the estimate of finite-size corrections to the dominant saddle point calculationcll. Actually, 
the evaluation of finite-size corrections in spin glasses (even in mean-field) is known to be quite difficult. Concerning 
the role of the cross-correlations between the noises we remind that the appearance of vanishing modes tends to be 
suppressed by cross-correlations JC < 1 which play the role (according to eq.(|lC])) of a thermal noise on the system 
described by the effective distance variables Zi = xi — yi. This explains why JC = 1 yields the lowest damage spreading 
temperature. We must stress that this transition should persist beyond mean-field. Actually the presence of a huge 
number of saddle points in phase space is not an exclusive mean-field feature but should persist in the presence of 
short-range corrections. Short-range corrections dramatically modify the height of the barriers but not their number 
which still could remain as large as in the mean-field approximation. 

To remark the advantage of considering Langevin dynamics for DS we have also investigated the DS in discrete 
glassy models without disorder such as the Bernasconi model with heat-bath dynamics. The difference between any 
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discrete dynamics and Langevin dynamics is that the type of correlator K in the last case may be chosen at will. So 
the case K, = 1 which has been studied for Langevin dynamics cannot be implemented for Monte Carlo dynamics. 
Each algorithm for Monte Carlo dynamics defines a given K. so the study of damage spreading in those cases remains 
more speculative because K. is essentially unknown. For the particular case of the heat-bath in the Bernasconi model 
we find a DS transition which separates two regions, a region T > T\ where is independent of the initial distance 
but finite and a region T < T± where D^, depends on the initial distance. T\ coincides (within numerical precision) 
with the mode-coupling transition Td- Although for heat-bath dynamics Doo depends on the initial distance we may 
not exclude that convergence to the asymptotic limit is extremely slow and convergence poor. 

We end up our discussion describing some opening problems. Much work on the DS has been devoted to the study 
of discrete dynamics in discrete systems the situation being the opposite for continuous-time dynamics. In this work 
we have shown that for Langevin dynamics we may choose specific noise correlations so it is easier to interpret what 
physical properties of the system we are looking at. It is difficult to adscribe any physical significance to properties 
measured for arbitrary discrete algorithms and this has been one of the major problems to interpret the large amount 
of numerical results obtained in the study of damage spreading in Monte Carlo simulations. Here we have seen that 
for Langevin dynamics such a task turns out to be easier. Still one of the major tasks which remain open is to 
understand better under which conditions there is a unique absorbing state for the damage. In another direction, one 
would like to get analytical proofs and equilibrium based analytical methods to investigate the connection between 
the DS transition and the topological properties of the phase space (such as the presence of saddle points in phase 
space). Finally, we would like to extend our research to real structural glasses where it could be interesting to study 
DS and investigate in which conditions such a transition is a precursor of the glass transition. 
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APPENDIX A: 



In this appendix we describe some of the main steps necessary to derive equations (|l9|)-(23). Consider two replicas 
{<r}, {t} submitted to the dynamics (pi]), 

^l = F i ({a})-na i + m , ^ = F i ({r})- m + r,' i (Al) 



where Fi is the force acting on the spin <x; due to the interaction with the rest of the spins eq. (|12[). TJae noises 
77,7/ are generically correlated < Tii{t)r]^{s) >= 2T/C({ct}, {r})S(t — s)Sij. Following Domany and Hinrichsenc3 we will 
suppose that the correlator K, is a generic function of the overlap due to the infinite-ranged nature of the model. So 
we will take /C({(t},{t}) a generic function of the equal-times overlap K,(Qd) whereO^ = jj^2iLi&iTi- Following 
the same steps as is usually done in the study of the dynamics of a single replica (seeEJ for details) we may write the 
generating functional for the dynamics in the Ito prescription, 

Z dyn = J[d<rdT}S(&i - Fi(a) + - ^(f,; - Fj(r) + nn - ry-) (A2) 

Introducing a new set of fields <7j,Tj and averaging over the noise we get, 

Zdyn — J [dadadrdr } exp(£(c, <r, r, f )) (A3) 

where 

C(a,a,T,t) = -i J dtY J vK°l-F l {a t ) + H ( r t i )-i J dt^M ~ F<(t') + /itf) - 

i ' i 

T J ^{^i-l+n f + ^i-l-rt) 2 ) (A4) 

i 

Because Zd yn = 1 we may average the dynamical partition function over the disorder. We use the cumulant 
expansion and retain only the first and second order terms exp(V) ~ exp (V + \{V 2 - (V) 2 )) where (.) stands for 
disorder average. Using eq. (ph we obtain the final result, 
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Zdyn = J [dadadTdr] exp(<S(er, a, r, f)) (A5) 

where 

S(<r,<7,r,f)=-i J dtJ2^m+^j)-i J dt^ijti + rfi-T J dt£(i±£(6* + f*) 2 + - f*) 2 ) 

2 2 2 

-~ | dt d S (^^|/'(C(t, S ))+^f|f//'(i)(i, S )) +^^f//'(Q(t, S ))+^f*a|/'(Q( S ,t))) 

i i i i 

I *^E(<^K^ /"(C(^s))+^f;r?rj s)) + ^f|r/ ^ /"(Q(i, *)) + f^ofrj f"{Q{s,t))) (A6) 

The correlation functions C,D,Q are defined by, 

c (*< s ) = ^ E < > < D ^ s ) = ^ E < ^ > - s ) = ^ E < ^ > ( A? ) 

z z i 

After introducing the necessary response functions, for instance R(t, s) = i < of cr| > we may solve the saddle 
point equations. Introducing also appropriate causality relations (for instance < ff<r| >= 0) we obtain the desired 
set of equations. Note that only the dynamical equation for the equal times overlap Qd{t) — Q{t,t) depends on the 
correlator /C. 
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In this appendix we prove the existence of the lower and upper bounds eq. 
To obtain the lower bound we start from eq. (23) and write, 



for T, 



o m 



1 dQ d (t) 

2 dt 



+ fj,(t)Q d (t) - TK(Q d ) - ^ / duR^^QP-^t.u)- 



p(p ~ 1) 



case a = 1, K = 1. 



duR(t,u)Q(t,u)C p - 2 {t,u) = 



(Bl) 



Let us assume now that we are in the high-temperature phase where the R and C are time-translational invariant 
and the fluctuation-dissipation theorem is satisfied, 



TR(t 



dC(t 



ds 



(B2) 



At first glance, this condition may seem too strong. Nevertheless, assuming FDT is allowed if we are in the high- 
temperature phase. At the DS temperature (which is deep in the high temperature phase, as we show in detail in 
sections III. A and III.B) the typical relaxation time for the quantities C and R is finite so we may introduce FDT 
with no harm. On the other hand, sect ion III.B shows that the DS transition is already present if we start from an 
equilibrium condition where condition (B2) is satisfied for all times. 

Using the inequality Q(t, s) < Qd(s) where Oswe obtain, 



\?^ + tit)Q*{t)-Ttc<g d )-£ 



* , dC p (t-u) ^ v _ u , p(p-l) f* , dC p (t-u)^ , , , 
du K - >-Q p >)- „ 1 / du y — J -Q d (u)<0 



du 



du 



(B3) 



Doing an integration by parts for the two integrals appearing in (B2) we may write, 

l^ir + Kt)Qd{t) T]C{Qd) ~ f (Q5 _1 w - er'co)^) - /* du dQP £ {u ' ) gp (* - u )) 



Pip - 1) 



(Qd(t)-Q d (0)C p (t) 



(i 



o 

t du dQM CP{t _ u)) < Q 



du 



(B4) 



Noting that Q d (t) is a monotonous increasing function of time (again, this is true only above the DS transition) we 
may write, 



t d u dQP t 1{u) CP(t-u)>0. 
du 



du 



(B5) 
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so, in the long-time limit t — > oo where the C(t) vanishes we get, 

l^ir + nQd(t) ~ Km) ~ f iQP ^ 1{t) ~ Qd{t)) - (B6) 

where we have replaced n(t) by its equilibrium value T + Finally, due to the monotonicity of Qd(t) we reach 
the desired inequality, 

< < T(K(Qd) - Q d (t)) + ^M^it) - i) (B7) 

If we approach the DS transition To from above we expect Qd(t) to relax very slowly to its asymptotic value 

Q d (oo) = 1. So, we expect \ < \ dQ g t (t) \ < (1 - Qd(t)) for large times. We may differentiate the inequality (|B7|) 

and put Qd — 1 which yields 



>i^ a ^ (r(1) _ 1)T+ fci>, (B8) 



At the DS transition the inequality is satisfied only if Tq > 2(i—K'{i)) ' ^ or case P ~ 2 the second member in the 



last expression of the inequality (BS) vanishes so the inequalities are never violated as soon as Qd < 1- This means, 
that there is no DS transition in the high-temperature phase for p — 2. The only possible transition occurs when 
C(t, s) and-4?(t, s) are no time-translational invariant anymore and this may happen only at Td. Actually, calculations 
of StariolcE3 show that the DS transition is present at the static transition temperature T s (which is equal to Td). 

As a curiosity (for which we have no analytical derivation) we note that this lower bound (for the case JC = 1) 
seems to coincide with the exact DS transition temperature for the fully asymmetric case a — (see section III.C). 

On the other hand, the upper bound can be obtained via a linear stability analysis of the equations ( pT| ) , ( |23| ) around 
the HT solution. Writing Q d (t) = 1 - e/(t) and Q(t,s) = C(t - s) - eg(t,s) with /(0) = 1 and g(t,t) = f(t) yields 
for equation ( p3| ) in the large t limit, 

i| = - mi -ru>>-f)/-*r^ ( -»>^ . (B9 , 



Finally, the inequality d9 g u u ^ > yields the upper bound. 
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